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Abstract 

Newton' viscosity law for the momentum flux and Fourier's law for the heat flux define Navier- 
Stokes hydrodynamics for a simple, one component fluid. There is ample evidence that a hydrody- 
namic description applies as well to a mesoscopic granular fluid with the same form for Newton's 
viscosity law. However, theory predicts a qualitative difference for Fourier's law with an additional 
contribution from density gradients even at uniform temperature. The reasons for the absence of 
such terms for normal fluids are indicated, and a related microscopic explanation for their existence 
in granular fluids is presented. 

PACS numbers: 



I. INTRODUCTION 



Granular fluids are of increasing interest to the Chemical Engineering, Physics, and Ma- 
terials Sciences communities, for different and complementary reasons. On the practical side 
are industrial issues of agricultural, pharmaceutical, and chemical significance for packing 
and transport of grains. More recently, the planned expeditions to the moon and Mars 
require an understanding of the surface regolith, a new form of granular matter. In spite of 
the growing phenomenology for granular fluids, the fundamental descriptions for the sim- 



plest states remain subject to question. An important c 



ass of questions involve the form 



Here, one new feature of the 



and conditions for a hydrodynamic description 
granular Navier-Stokes hydrodynamic representation is addressed. 

There are two peculiarities of the Navier-Stokes equations for granular fluids relative 
to those for normal fluids. One is a dissipative source of energy due to the inelasticity of 
collisions among the grains, easiliy understood as a consequence of the mesocopic nature 
of the constituent particles. More puzzling is the modification of Fourier's law, with an 
additional contribution to the heat flux s due to a density gradient 

s = -\VT - fxVn. (1) 

Here A is the thermal conductivity and fi is the new transport coefficient characterizing heat 
flow in an isothermal, inhomogeneous flnid. Althonglr plrysical interpretations ot this latter 

effect have been given [6], it is instructive to associate it with the fundamental differences 
between normal and granular fluids. At the continuum level fluid symmetry alone leads 
to the general form ([1]) in both cases. Therefore, to better understand the presence of the 
density gradient contribution for granular fluids, it is useful to ask why it is absent for normal 
fluids. 

This question is answered here at the level of the Boltzmann kinetic equation for a gas 
and using the more general statistical mechanics of linear response. In fact, a definitive 
conclusion requires information beyond continuum mechanics from the microscopic basis 
for hydrodynamics. The hydrodynamic fields are averages of the densities associated with 
the global invariants of mass, energy, and momentum. The hydrodynamic equations are 
then obtained from representations of the corresponding fluxes in the average microscopic 
conservation laws for these densities. There are two central observations here that are 
responsible for /x = in normal fluids. The first is that the mass flux is equal to the 



momentum density. The second is that the reference state (equihbrium or local equilibrium) 
is entirely determined by the invariants. The first implies that there is no dissipation in the 
continuity equation for the mass density. The second imposes important constraints on 
the representation for transport coefficients in terms of correlation functions for the fluxes. 
The most compact version of these is Onsager's theorem |^], which is discussed next. Then 
more explicit realizations of these correlation functions are obtained from low density kinetic 
theory and formally exact statistical mechanics. 

It should be emphasized that the existence of the coefficient /i is not in question, only its 
origin and interpretation. The coefficient has been determined for idealized isolated cooling 
granular gases ^, and measured in both simulations 9| and experiments |10j] for shaken 
fluids in a gravitational fleld 

It is a pleasure to dedicate this work to Keith Gubbins whose contributions to the kinetic 
theory and statistical mechanics of fluids have influenced me for more than thirty years. He 
has been an exceptional combination of mentor, colleague, role model, and friend to me and 
many others. 



II. A CONSEQUENCE OF ONSAGER'S THEOREM 

As noted above, the Navier-Stokes conditions of fluid symmetry and small spatial gra- 
dients do not constrain the "constitutive equation" for the heat flux beyond the form ([T]). 
However, in 1931 Onsager made a seminal observation relating the transport coefficients of 

Q 

a normal fluid to the underlying statistical mechanics of the fluid at equilibrium [7|. The 
Onsager regression hypothesis states that on long space and time scales the decay of sponta- 
neous fluctuations in an equilibrium fluid is governed by the same laws as for nonequilibrium 
states displaced slightly from equilibrium. This is effectively what has become formalized 
as linear response theory. The symmetry property of these equilibrium fluctuations result- 
ing from microscopic time reversal invariance of the dynamics leads to relations among the 
transport coefficients. These are the familiar Onsager relations for a fluid mixture. However, 



it is less we 



in Eq. (HD Ul, ll2| 



1 recognized that these same relations for a one component fluid imply that n = 



To show this explicitly, consider the exact macroscopic balance equations for the number 
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density n{r,t), energy density e(r, t), and the momentum density g(r, t) 

dtn{r,t) + m^^V ■ g(r,t) = 0, 



(2) 



dte{r,t) + V ■s{r,t) =w{r,t), (3) 

dtgiir,t) + djUj{r,t)=0, (4) 

where g(r, t), s(r, t), and tij{r,t) are the associated mass, energy, and momentum fluxes 
and m is the mass. Here, it has aheady been assumed that the flux of mass is the same as 
the momentum density. In fact, this requires proof from the underlying microscopic balance 
equations. Also, w{r,t) is an energy source term that could be due to an external force 
doing work on the fluid, or the internal collisional energy loss of a granular fluid. The fluxes 
have a contribution due to convection and a purely dissipative contribution that occurs in 
the local rest frame at each point of the fluid. This separation is easily identifled by a local 



Galilean transformation with the results 11| 



g(r,t) = mn(r,t)u(r,t), (5) 

Si{r, t) = (e(r, t) + p(r, t)) M,(r, t) + t*^ (r, t) Uj{r, t) + s*(r, t) (6) 

Uj (r, t) = p{r, t)6ij + mn(r, t)ui{r, t)uj{r, t) + t*j (r, t) . (7) 

Equation ([5]) deflnes the local flow velocity u(r,t), while ([6]) and ([7]) deflne the irreversible 
energy flux s*(r, t) and momentum flux t*-{r,t) in the rest frame. Finally, p{r,t) is the 
hydrostatic pressure which must be specifled as a function of the local energy and density. 

The set of exact equations ©-(Ej) become a closed set of hydrodynamic equations once 
the "constitutive equations" for s*(r, t) and t*j (r, t) are given in terms of the number, energy, 
and momentum densities. In practice, it is useful to introduce a conjugate set of variables 
through the change of variables 

Xair, t) ^ -^il^l^, y„(r, t) ^ (n(r, t), e(r, t), g(r, t)) . (8) 

Here s ({^/^(r, t)}) is the equilibrium entropy density for the fluid as a function of the number, 
energy, and momentum densities at each point of the fluid. For example, it follows that 
X2 = —l/T which deflnes the temperature as a function of the density and energy. Note 
that although this entropy function is deflned for an equilibrium fluid, it is used here simply 
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as the mathematical generator of a change of variables. As such, it applies even to states 
far from equilibrium and also to the granular fluid. 

For states near uniform equilibrium, the leading contributions to s*(r, t) and t*j (r, t) are 
linear in gradients of the {i/a}, or equivalently, gradients of The resulting equations 

are the Navier-Stokes order hydrodynamics 

dtya{r,t) + V ■ ja{r,t) = Wa, (9) 

where the fluxes are given by 

j«(r,t) =}^^Hr,t) +j:(r,t), j,(r,t) ^ (g(r, t), s(r, t), t,,(r, t)) . (10) 

Here, ja°^ are the Euler order contributions identified from ([S]) and ([7]), and j* are the irre- 
versible contributions. At Navier-Stokes order the latter are given by the linear constitutive 



equations 
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j:(r,^)=-EW({l//3(r,t)}) : (11) 

The energy flux in this approximation is the generalized Fourier's law ([T]) and the momentum 
flux is the usual Newton's viscosity law. 

Up to this point only the macroscopic conservation laws, gradient expansion near the ref- 
erence homogeneous state, and fluid symmetry have been used. For a normal fluid Onsager's 
regression hypothesis allows identification of the transport coefficients Laj3 in the form 

POD 

^a/3= / dtC^f,{t), (12) 

where Caisit) is an equilibrium time correlation function for two microscopic fluxes corre- 
sponding to the macroscopic j^. These fluxes have a definite parity under the transformation 
t — i> —t, and reversal of all particle velocities {vj}-^{— Vj}. Since both the dynamics and 
the equilibrium ensemble are invariant under this transformation it follows that 

Ca/sit) = TaTfsCpait), (13) 

where = ±1, depending on the parity of the associated flux. Therefore Lq,^ has the 
symmetry 

La/3 = TaT/sL/Sa- (14) 

This is Onsager's theorem. 



The coefficient in ([T]) is identified as 

fi = L2i-^\T. (15) 

An important observation now is that the macroscopic conservation law for the number 
density, the continuity equation, has no dissipative contribution, Li^ = 0, since the mass 
flux is equal to the momentum density. Consequently, Onsager's theorem gives 

V = Li^ = 0, (16) 

and Fourier's law becomes 

s*(r, t)= -X {{yp{r, t)}) VT(r, t), A = L22T-2. (17) 

(There is no coupling to the velocity gradients, Las = 0, from fluid symmetry). Thus, /i = 
for a normal fluid as a consequence of the fact that the number flux in the continuity equation 
has no dissipative contribution, and the symmetry of the correlation functions resulting 
from the fact that the reference state is a function of the invariants. These properties are 
demonstrated more explicitly in the further microscopic elaboration below. 

The above analysis holds as well for granular fluids, except for Onsager's theorem. The 
condition L1/3 = still applies but the key symmetry, (fT3l) . rests on both the equilibrium 
Gibbs reference state and the invariance of the dynamics for the system. Neither the ho- 
mogeneous reference state nor the dynamics of a granular fluid have these symmetries, so 
neither the usual form of Onsager's theorem nor the conclusion that fi = can be extended 
to granular fluids. This failure of Onsager's theorem occurs as well for normal fluids in 
nonequilibrium stationary states. However, it is possible that some other symmetry could 
apply to enforce the usual form of Fourier's law for granular fluids. Thus, it is important 
to look in more detail at the mesoscopic (kinetic theory) and microscopic (statistical me- 
chanics) representations of this constitutive relation. This is the objective of the next two 
sections. 

III. KINETIC THEORY 

The simplest fluid is a low density gas for which the appropriate kinetic theory is given 
by the Boltzmann equation. The usual derivations of this equation for a normal gas can 

6 



)e extended to the granular gas as well, to account for binary collisions that are inelastic 
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16| . It provides an important testing ground for the derivation of hydrodynamics and 



for exploration of conceptual issues as well. The application to a normal gas is considered 
first, and then its extension to a granular gas. 



A. Normal gas 

The Boltzmann equation for the one particle reduced distribution function /(r, v,t) is 

(5, + vV)/ = C[/,/]. (18) 



where C[f, f] describes uncorrelated elastic binary collisions llj. The notation denotes that 
C[f, f] is a bilinear functional of /. Its detailed form will not be required here, only the fact 
that particle number, energy, and momentum are conserved as represented by the properties 

J dwa^ClfJ] = 0, a«(v) ^ (^l, imt;^mv^ . (19) 

The set of functions are known as the summational invariants. Their averages are also 
the hydrodynamic fields of the last section 

y^ir,t) = j rfva,(v)/(r,v,t). (20) 

The macroscopic balance equations (I2l)-(j4l), or in the equivalent compact form ([9]), follow 
from this definition by differentiation with respect to time and application of the Boltzmann 
equation. In this way the fluxes are identified as 

ja(r,t) = y civb„(v)/(r,v,t), b«(v) = va„ (v) . (21) 

Hydrodynamic equations result from these exact consequences of Boltzmann's equation 
when the solution / approaches a "normal" form on some length and time scale, expected 
to be long compared to the mean free space and time 5|. A normal distribution is one for 
which all space and time dependence occurs through the hydrodynamic fields 

/(r,v,t)^/(v|{y4). (22) 

The notation /(v | {ya}) indicates a functional of the fields ya{r,t) throughout the system 
(equivalently, and for practical purposes, it is a function of the fields and all their derivatives 
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at the point of interest). An example of a normal distribution is the local equilibrium 
distribution 

hiV) ^ n (^) exp {-mV'/2T) . (23) 

where V = v — u(r, t) and in all of the following, units are used such that Boltzmann's 
constant ks = 1- This distribution is parameterized by five fields n, T, and u which, as the 
notation suggests, are chosen to be the same as the hydrodynamic fields defined in fl20|) and 
([8]). This is expressed by the condition 

J dwa^ if - fe) = 0. (24) 

The local equilibrium distribution function depends on the fields, but not their gradients. In 
fact, it is a solution to the Boltzmann equation to zeroth order in the gradients, as follows 
from the second important property of the collision operator 

C[fe,M=0. (25) 

More generally, the normal solution to the Boltzmann equation can be constructed as an 
expansion in the gradients with as the leading order contribution. This is done by the 



familiar Chapman-Enskog procedure and is carried out in the Appendix. 

It is useful to provide a geometrical representation of the normal solution. First, define 
a set of functions "conjugate" to the {aa} by 

iljy = -— , / dvaaipu = 5au- (26) 



The second equality follows from fl20|) and fl24l) and shows the sense in which {aa} and {ipa} 
form a biorthogonal set. Next, define the projection operator whose action on an arbitrary 
distribution function h is 

Vh = ip^ / dva^h. (27) 



It is easily verified that V has the property of a projection operator = V. The utility 
of these definitions is the decomposition of / into its local equilibrium distribution plus a 
remainder that is in the orthogonal subspace defined hjQ = l — V 

f = fi + 5f = fi + V6f+ Q6f = fe+ Q6f. (28) 

The last equality follows from (^^, V6f = 0. Since fi is a function of the fields, and not 
their gradients, all contributions to a normal solution due to gradients must come from 
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Sf = QSf. Furthermore, since fe is a solution to the Boltzmann equation at zeroth order in 
the gradients (see fl25|) ). 6f = Q6f is at least of first order. This allows a decomposition of 
the fiuxes (121 p into Euler and irreversible contributions as in the previous section with the 
identifications 

jW(r,t) = J dvhje, j:(r,t) = J dvKQ5f. (29) 

An immediate consequence of the property bi (V) = as (V) is = 0, resulting in the 
continuity equation for the number density. 

The Boltzmann equation determines the detailed form for 6f. Assuming a normal form, 
the Boltzmann equation to first order in the gradients becomes (see Appendix) 

M/ = -(Q7a)-VrXa, (30) 

where L is the linear Boltzmann collision operator 

Lh^-C[fe,h]-C[hJe]. (31) 
A transformation to the conjugate variables has been made using 

VrZ/a = Oau'^rXu, Qua = 9au = (32) 

and 7a (v) are the conjugate fiuxes are defined by 

7„(v) = Y^j^ (v) g^^ = v^l^ (33) 

The appearence of the orthogonal projection Q on the right side of fl30l) assures that so- 
lutions to this equation exist. This is given by the Fredholm alternative for such linear 
inhomogeneous equations which states that the right side must be orthogonal to the null 
space for the adjoint of C, in this case given by VC = 0. The formal solution can be written 

poo 

Sf = -Q t^te-^*Q7a(V)-VrXa, (34) 
Jo 

In general, an arbitrary solution to the homogeneous equation L6f = could be added to 
this, but the condition (IMl) . V6f = excludes such terms. 

Use of this formal solution in (!29l) gives the linear constitutive equations (fTTj) . and the 
transport matrix Lafs is identified as 



oo 



drCaf^ir), (35) 







Co.(3{r) = j dw (Qtb,) e-^^Q7^. (36) 

where Q} is the adjoint of Q. This makes exphcit the low density form for the general 
phenomenological postulate. It is given in the form of a Green-Kubo expression, where the 
"time correlation function" CQ,/3(r) is a flux - conjugate flux correlation function. Further 
simplifications are possible by noting that both the conjugate densities ip^ and conjugate 
fluxes can be written as linear cominations of the densities and bj^, respectively 

'ipu = fi"^ Cuada, 'lu = fl"^ C^^^K- (37) 

a a 

Note that the (adjoint) projection operator implies that Ci/3(r) = since bi oc as (the 
mass flux is the momentum density) and projects orthogonal to the set of densities \aa\- 
The [I coefficient vanishes for similar reasons. It is given by 

/x= r drC^^ir)^-^ \t . (38) 



The conjugate flux 71 is proportional to the conjugate density 

71 = V- — |t= ^ It W-i- (39) 

Since Q projects orthogonal to the conjugate densities {V'a}, i-e. Qt/'s = 0, the projected 
flux Q7/3, the correlation function C2i(t), and consequently /i all vanish 

T dxi 



^ 1^ /"^^ / (2^b„)2e-^^QV'3 = 0, (40) 



mn dn 

In summary, the correlation functions determining the transport coefficients involve fluxes 
projected orthogonal to associated densities. The fluxes associated with the density are 
themselves densities, and hence such correlation functions vanish. In particular those asso- 
ciated with dissipation in the continuity equation and the contribution to the heat flux from 
density gradients are zero. 



B. Granular gas 

The analysis for the granular Boltzmann equation proceeds in a similar way [XT'], and the 
details are given in Appendix A. The correlation functions for the transport matrix L^^ in 
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this more general case are 

Ca^(r) = jdY (Qtb,) e-^- (e^^Q7)^ • (41) 

C^l^V^-V + L, e = (42) 

The parameter ^ is called the cooling rate since it determines the relative rate of change 
of the temperature in the HCS due to inelastic collisions. The densities {oq} and fluxes 
{bo,} are the same as those in (fT9l) and (12T!) . The projection operator is given by (12T|) but 
with new conjugate densities {ipa} and fluxes {70} defined in terms of the reference local 
homogeneous cooling state (HCS) f^^^ instead of the local equilibrium state 

V'. = 7« = (43) 

The reference HCS is the normal solution to the Boltzmann equation in the absence of 
gradients 

^^Vv(V/W)=C[/(°),/(°)]. (44) 

The cooling rate ^ vanishes in the elastic limit and f^^^ — > fe in this limit. 

Aside from these differences in the reference state for granular fluids, the generator for 
dynamics in (jHj) shows significant differences. The linearized operator L is now that asso- 
ciated with the Boltzmann collision operator for inelastic collisions. In addition, there is a 
velocity scaling operator ^^Vv • V (the derivative operates on everything to its right) which 
compensates for the cooling generated by the collisions, as illustrated in the exact balance of 
these effects in (144|) . Finally, there is the matrix A whose explicit form is given in ( ]A22[) and 
(IA24I) of the Appendix. The relevant point here is that its eigenvalues (O, — |^) are the 
same as the smallest eigenvalues of C. Thus, the entire generator for the dynamics has a null 
space, just as for normal fluids, and the projection Q assures that the fluxes are orthogonal 
to the invariants defining this null space. This is a necessary condition for the existence of 
the integral defining the transport matrix Laf3- 

It is seen that the condition for the continuity equation still holds Ci^(r) = 0, for the 
same reason as in normal fluids, Q^bi oc Q^sl^ = 0. However, the related conditions for the 
transport coefficient fi to vanish are no longer satisfied in general. The relevant correlation 
function C2i('r) is no longer determined by the single flux 71 but is coupled to 72 as well. 
Neither of these fluxes in simply proportional to the conjugate densities {ipa} and hence the 
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action of Q on each is non zero. This difference occurs for granular fluids because f'^^^ is no 
longer determined entirely in terms of the {aa} and hence there is no simple relationship 
between the sets {aQ:},{ba} and {V'ojjlTa} as expressed in ( 1371) . The explicit form for /x is 
obtained in Appendix A 



dn 



_2|l|,(e---..-e-.-.V)s(v^|^|„)} 



(45) 



The new conjugate densities {V'a} are invariants of the new dynamics generated by {C — Aa). 
The projection operators again assure that there is no contribution from the invariants 
which is a necessary condition for convergence of the r integral. In general, however, neither 
vdf^^ydn \t nor vdf^^yde |ris a linear combination of these invariants and the action of Q 
on them does not vanish. Only in the elastic limit does Qvdf^^ydn It become proportional 
to QVa = 0- III this limit A2 = Ai = and the coefficient of \df^^^/de \t also vanishes, 
confirming /i = for a normal fluid. 

IV. FORMAL LINEAR RESPONSE 

The analysis of hydrodynamics from kinetic theory can be generalized by the formal ap- 



plication of nonequi 
in references 
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ibrium statistical mechanics to granular fluids. The details are described 



19| . The starting point is the Liouville equation for the particle phase 



space density p (T, t), where T = {qi, ..q^v, vi, ..vjv} denotes a point in the 6A^ dimensional 



phase space 20|]. First, a homogeneous normal solution to the Liouville equation is ph (r) 



identified, representing the homogeneous cooling state (HCS) for an isolated system 



- 1 ^ 

Cph = 0, CX = -a Yl ■ (V,X) + LX. (46) 



1=1 



Here L is the Liouville operator for A^ hard inelastic spheres and (h is the associated cooling 
rate in the HCS. In addition, there is a scaling operator |C/iVv, ■ Vj for each particle. 
Clearly, C is the A^ particle generalization of the kinetic theory generator C of fH2|) . Next, 
small spatial perturbations of this state are induced through an associated local HCS peh- 
The response of the hydrodynamic fields at a later time due to these initial perturbations is 
characterized by response functions, represented as time correlation functions composed of 
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the particle phase functions for the hydrodynamic fields (corresponding to the of the 
last section) and the functional derivatives of pih with respect to the conjugate fields {Xa} 
(corresponding to the ipa of the last section). Since these response functions must exhibit 
hydrodynamic excitations at long times and long wavelengths, the transport matrix Lq^ can 
be identified. It has the representation in terms of time correlation functions again, as in 
(135|) but now the results are formally exact without the restrictions of the kinetic theory. 

For the purposes here it is sufficient to display only the final result for the coefficient /x. 
It has a form similar to that of (l38l) 



/i = /io + ^"rfrC2i(r)^ It. (47) 

There is an additional term /ig that does not have the form of a time integral of a correlation 
function, and is due to both the singular hard sphere dynamics and the dissipation. It 
vanishes in the low density limit, and so does not appear at the level of the Boltzmann 
equation. It is given by 

^'o = ^v-^ j cirs(r)-M(r). (48) 

where S(r) is the volume integrated phase function for the energy flux (corresponding to b2 
in the kinetic theory analysis), and M(r) is the space moment for the functional derivative 
of pih with respect to density 



M=/*--(i^t) . (49) 

5y=0 



6n (r) 



It is easily seen that fiQ —>■ in the elastic limit, for which pih becomes the corresponding 
equilibrium ensemble. Otherwise it is non-zero. The correlation function in P7|) has a form 
similar to that of (HTl) 



C2i(r)^ |t= V-'I dV (Qts) e-^^ {e^^QT)^ (50) 
where the conjugate flux is 

T„ = -(£-A)M, (51) 

Here Mq, is the space moment of the functional derivative of pih with respect to ya holding 
the other {y/^} constant, and evaluated at the HCS. 

The projection operator Q projects orthogonal to the invariants of the dynamics gener- 
ated by (£ — a) . Consider now the elastic limit for which pih becomes a local equilibrium 
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ensemble pie- To be specific, let that be the local grand canonical ensemble. The correlation 
function simplifies to 

C2i(r)^|T -V-^ j dV (Qts) e-^^QlM, (52) 

M ^ -pem j drr j dr' n (r') , (53) 

where u is the activity, pe is the strict equilibrium ensemble, and n (r) is the phase function 
corresponding to the number density. Then using Ln (r) = — m^^ V • g (r) , where g (r) is the 
phase function representing the momentum density, the flux LM is found to be proportional 
to the total momentum P 

IM = p,m [ drr [ dr' ( ^JlMu^ Ln(r') 
J J \6n{r) Jsy^^ 

Consequently, QLNL = and C2i(r) = 0. It was already noted that po vanishes in the 
elastic limit, so /i = as well and the usual form of Fourier's law is recovered. 

More generally, for inelastic collisions po (except at low density) and none of the 
fluxes T are simply proportional to the invariants. Consequently, QT ^ and p ^ 0. 



V. DISCUSSION 

The theoretical "discovery" that Fourier's law for a granular fluid has an additional 
term proportional to the density gradient provides an interesting qualitative difference from 
normal fluids. There has been much discussion about this additional term and attempts to 
detect it in simulations or experiments. The theoretical analysis here provides a different 
perspective, in hindsight, that the surprising nature of this difference between Navier-Stokes 
hydrodynamics for normal and granular fluids is the absence of this term in the former case 
rather than its presence for the latter case. This circumstance is similar to the discovery of 



generic long range correlations in nonequilibrium states, absent at equilibrium [2l|]. Both 
provide examples of the very special balance of competing effects for the equilibrium state, 
in contrast to the qualitatively different behavior for generic nonequilibrium states. 



14 



The vanishing of the transport coefficient for the contribution to the heat flux from a 
density gradient follows from the simple structure of the local conservation law for the 
number density, and the characterization of the equilibrium state in terms of the dynamical 
invariants. The number conservation law relates the density to its flux. However, this 
flux is itself proportional to a density for one of the invariants - the momentum. The 
transport coefficients are time integrals of correlation functions composed of these fluxes and 
corresponding "conjugate" fluxes generated from the equilibrium ensemble. Convergence of 
the time integrals requires that there be no time independent parts to these correlation 
functions. Thus, the correlation functions are constructed from those parts of the fluxes 
that have their invariant parts subtracted (or projected) out. As the number flux is itself 
an invariant all transport processes coupling to the number flux will therefore have zero 
transport coefficients. This leads directly to the vanishing of all dissipative contributions 
to the continuity equation. In addition, since the equilibrium ensemble is a function of the 
invariants it results that the conjugate fluxes are linear combinations of the fluxes in the 
conservation laws. Thus any contributions from the number flux, itself a density, has no 
remainder once its invariant parts are subtracted out. This leads to the vanishing of the 
coefficient fi in Fourier's law. 

The simple structure of the conservation law for number density does not depend on the 
state of the system considered. Thus, there are no dissipative contributions to the number 
flux for any state, equilibrium or nonequilibrium. In contrast, the conjugate fluxes depend 
in detail on the reference state about which hydrodynamic excitations are being considered. 
The simple relationship of the conjugate fluxes to the fluxes in the conservation laws for 
equilibrium states cannot be expected more generally for any nonequilibrium reference state. 
This is the case for both normal and granular fluids. An example of the former is the 



hydrodynamic excitations about uniform shear flow 13| where Fourier's law in the form of 
Eq. (P) applies with 7^ 0. 

In closing, it may be useful to display concrete expressions for A and yU in a granular gas 



obtained 
equation 



rom an approximate evaluation of the expression given above from the Boltzmann 



221 for ci = 3 



, nT 5 T ^ 2i 

A = — — — — -, /i = — A- 



m 2(z/-20' n 2Z/-3C 

where ^ = Sz/q (1 — ct^) /12 is the cooling rate of (jHj) calculated in this same approximation, 
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uq (1 + a) (49 — 33a)/48, and vq = IGna'^ (7rT/m)^^^ /5 is an average collision rate. In 



the nearly e 



astic limit A IbnT/Amvo and {1 — a) 75T /16mz/o. The simulations of 



references js] and 8|] confirm the more general theoretical prediction over a wide range of 
values for the restitution coefficient a. 
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APPENDIX A: NORMAL SOLUTION TO BOLTZMANN EQUATION 

In this Appendix the normal solution to the Boltzmann equation for a granular gas is 
obtained up through first order in the gradients jl?! . 

/ = /(°)+F,-Vrl/a + " (Al) 

The leading order term f^'^^ and coefficient Fq are functions of the actual hydrodynamic 
fields {Ua}- As a normal solution obeys a condition analogous to (l24l) 

J dva^ (/ - /(°)) = 0. (A2) 

Only the case of hard sphere interactions is considered, so there is no internal energy scale. 
Then, from dimensional analysis, they have the forms 

/'°'-«(^)''"/'°'-(n, , V=y^(v-U), (A3) 

F.^{^fFHy), F,^-1(-)"'^f:(v). F3 4(i^)^F;(V), (A4) 

The temperature T for hard spheres is related to the energy by e = |nT + ^mnll'^, and the 
asterisk denotes a dimensionless quantity. Substitution of (I A II) into the Boltzmann equation 
gives 

^ Vr2/„ + • V + - (%. + v-VrZ/.) = C[f, /], (A5) 



and using the macroscopic balance equations ([9]) the time derivative can be expressed in 
terms of the gradients 

QfiO) QY \ 

'iiT^'df' + ■ V + ■■ j (ti;. - V ■ + vVrZ/.) = /]. (A6) 
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To zeroth order in the gradients this equation determines /^"^ 

^-^w^[f'\f'^]=C[f'\f'\ (A7) 

where the bihnear functional dependence of the energy loss w^j has been made explicit to 
show that here it is evaluated to lowest order. Since the source Wa occurs only in the energy 
equation the left side of this equation can be made more explicit 

[/'»'. /'°'] = [/'"'] = [/'"'. /'"'] = i« (<i + V . V) /<»), (A8) 

and the cooling rate ^ = —2w \_f^^\ /*-°-'] /3nT has been introduced. Equation (jATj) for /'^^^ 



becomes 



^eVv(V/W)=C[/W,/W]. (A9) 



Although simple analytic forms for the solution to this equation have not yet been found, 
its behavior for small and large velocities is known, and good approximations exist more 
generally. Furthermore, it has been studied numerically using Direct Simulation Monte 
Carlo. For our purposes, therefore, it can be considered as known. In the elastic limit the 
solution to C[/(°), /W] = is the Maxwellian (|23l). 

To first order in the gradients Eq. (1A6P determines Fq, as the solution to 



, 9Fa dw 
LF^ + ( w—- + — F2 

oe dya 



VrVa = (V- Vrl/. - Vr " ji"^ (AlO) 



Lh = -C[fe,h]-C[h,fe]. (All) 

It is understood here that w = Wa [f'^^K f^^^ ■ The Euler order flux is defined in terms 
of the lowest order distribution as in (12^ 

jW(r,t) = J rfvb„(v)/(°)(r,v,t). (A12) 

It follows directly that the second term on the right side of (lAlOp can be written 

^ Vr ■ (r, t) = ^^j rfva, (v) (v) ■ V^Va = V (v^, (v)) ■ V^Va, (A13) 

where ipa and the projection operator V are defined mutatis mutandis as in ( l26i) and ( 127|) 

5/(0) 



j dvtta (v) (v) = Sau, Vgiv) = tpuM j dwa^ (v) g (v) . (A14) 
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The gradients of Ua are arbitrary in (lAlOp so their coefficients give the desired equations for 



w— + £ j F„ + = -Q (vV^. (v)) , Q = l-V. 



(A15) 



This can be simphfied further by noting that w oc {nTY^'^n^^'^ and using the scahng of 
(lAl 

d 2w d dw Iw dw 2 dw w ,,^„s 

(A16) 



dw 1 w 



w 



to get 



de dn dT ' On 2 n 



3r, 



dw 2 9w w 
de 3n dT nT ' 



The constants Aq, 



(£ - A,) - 6^1^— F2 = -Q (v^„ (v)) 
£ = i^Vv V + L. 



are eigenvalues of the operator ^^Vv ■ V + £ 



(£ - A„) <~Pa = 



(A17) 
(A18) 

(A19) 
(A20) 



The eigenfunctions (fa are hnear combinations of the set {V'a}- Since the operator Q projects 
orthogonal to this null space the Fredholm alternative is satisfied for these integral equations 
and their solutions exist. 

The set of equations (1A17P can be written in matrix form 



{CI - A) = -Q (v^, (v)) = -g-'^Q (7/3 (v)) , 



with 



w 



— Aa^a/J — — 5al5/32; 



(A21) 



(A22) 



where / is the identity matrix. Also, the fluxes 7/3 are defined as in ([33] 

9/(0) (v) 



7a(v) = V(V^C/)^ 



dXc 



(A23) 



Next, perform a simlarity transformation with the symmetric matrix g to get the form 



{CI - A)„^ {gF)^ = -Q{la (v)) , A = gAg-\ 



The solution flAip can be given the representation 



Sf = f-f 



- f(o) = _ 



/■oo 
. ^0 



(A24) 



(A25) 
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The irreversible fluxes are obtained from (129 



The transport matrix is given by (!36|) with the correlation functions [23 1 

Gap{t) = Jdv (Q^K) e-"^ {e^'Ql)p • 
In particular, the correlation function determining the transport coefficient is 



C2i(t) = y rfv (Qtba) e-"^ [e'^'Ql), 



A more explicit form is obtained by direct solution to the equations for Fi and F2 



^1 - ^T7^ rT^2 



■4(A2- Ai; 

¥2 = - dre-(^-^^)"Q(v^2(v)). 
Jo 

The coefficient /i is then 

/i = drC,,ir)^ \t= J rfv (Qtb2) (f, + F2|^ |t 
Substituting (lA29|) and (lA30D leads after some rearrangement to 



(A26) 
(A27) 
(A28) 

2 

(A29) 
(A30) 

(A31) 



/i = dr j dv{Q%2) |e-(^-^i)"Qv 

)q(v 



(£-A2)t _ g-(£-Ai)r 



9e 



(A32) 
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